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Abstract 

The dynamics of one and two pointlike vortices in a planar quantum gas of spin-0 particles 
confined in an annular region is considered. New analytical and numerical solutions are found. The 
concept of stationarity radius, related to the doubly connected nature of the annulus, is defined. 
It is seen that the existence of these radii has great impact on the behaviour of the vortices. It 
is shown that, because of the existence of the stationarity radii, vortices exhibit similar behaviour 
regardless of the sign of their winding number. The energetically stable vortex solutions are studied 
qualitatively. 
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I. INTRODUCTION 



Several systems in low temperature physics exhibit different kinds of vortex solutions. 
Because these solutions have many interestingproperties, vortices have been extensively 
studied both theoretically and experimentally . In the current paper the 

starting point is a quantum system of spin-0 particles in the plane. This quantum system 
could in reality be, for example, a flat ordinary BEC or superfluid "^He-film. These systems 
can be described by an order parameter field taking complex values. The system under 
study has a global f/(l)-symmetry and thus the first homotopy group, classifying pointlike 
topological defects is isomorphic to the additive group of integers. The vortices are pointlike 
objects characterized by a single integer, called the winding number or topological charge. 

The pointlike nature of the vortices described above can be described as point particles in 
a well known way (for example jsl,!^). We will use this particle-description. As a consequence 
the effects of the vortex core are left out of the discussion. This is an approximation that 
can be justified when the vortex size of the vortex core is small compared to the dimensions 
of the whole physical system and the distances between vortices. This is in many actual 
cases true, and thus describing planar vortices as point particles presents a rather good 
approximation of real planar vortices. 

Our main interest is in a homogeneous system confined in an annular region (ring). We 
will incorporate the appropriate boundary conditions by using a simple Green's function 



method arising from the hydrodynamical description of the gas (eg. [UJ, ll2||). The choice 
of an annular region as the object of study is motivated by two reasons. First, the annulus 
is the most simple doubly connected domain in and the Green's function needed in our 
analysis is quite easily expressed analytically {l^. This enables simple numerics and also 
some analytic observations of vortex behaviour inside the annulus. Also, using the theory of 
conformal mappings, the annulus can be mapped to other doubly connected regions in M^. 
The second reason, more practical in nature, is that a ring shaped BEC can be produced 
by creating the condensate in a trap characterised by the so called "mexican hat" -potential. 
This means that experimental studies of vortex behaviour in an annular BEC are possible. 

The dynamics of point vortices in an annulus has previously been studied in [l3|. We 
obtain the same results along with new results. We will show that there always exists 
stationary (though not energetically stable) one- and two-vortex solutions in the annulus. 
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These solutions are characterized by specific values of r (the distance from the origin in polar 
coordinates) which we will call the stationarity radii. These radii divide the annulus into 
two subsets in the radial direction in such a way that position of the vortices with respect 
to the stationarity radii greatly affects the characteristics of a particular vortex solution. 
This fact has, at least to our knowledge, been unnoted before and is responsible for the new 



13|. The stationarity radii can, in general, be seen as 



solutions that are not mentioned in 

relatives of the equilibrium points (or critical points) of the Green's function in the annulus 



studied in eg. 



161, Q. They are however not the same thing, because the calculation of the 
stationarity radii takes into account the self-energies of the vortices. These self-energies are 
associated with the harmonic part of the two-dimensional Green's function. 

Demonstrating the importance of the stationarity radii, we find analytic solutions describ- 
ing corotating (the angular velocities of vortices have the same sign) and counterrotating 
(the angular velocities of vortices have opposite signs) motion for both vortex-vortex and 
vortex-antivortex cases. These solutions could not exist without the existence of the station- 
arity radii. With numerical simulations we further demonstrate how the stationarity radii 
along with the initial positions of the vortices play an essential role in the characteristics of 
the solutions of vortex motion. 



II. VORTICES AS ELEMENTARY PARTICLES IN A PLANAR QUANTUM SYS- 
TEM 

A. Description of planar vortices as point particles 

We will briefly review a simple way by which the point-particle theory of planar vortices 
(e.g. j^, y, [l^) can be derived by separating the vortices from the non-vorticical fluid part. 
The starting point is the semiclassical Lagrangian for a two-dimensional superfluid consisting 
of bosonic, spinless particles 

L[0,0*] = j d2a;Q(0*9o0-9o0*0)-^|V0p 

-f/(a3)0*0-V(0*0)). (1) 

The function U represents external potentials and V describes particle-particle interactions. 
Instead of the fields and 0*, the gas can also be described by hydrodynamical variables p 
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and V which are the particle density of the fluid and the velocity field of the particles. 

We will next choose the following parametrization, a two-dimensional analogue of the 
usual Helmholtz-decomposition of three-dimensional vector fields, for the velocity field v: 

v' = 9V + e'^djA. (2) 

Function if corresponds to the regular part of the phase of while the function A is associated 
with the irregular (singular) part of the phase. The familiar quantization condition for 
pointlike vortices with winding numbers {q^, ■ ■ ■ q^} in the system leads to a simple 
expression 

N 

Aix) = -Y,<l''Gix,x^), (3) 



a=l 

where the two-dimensional Green's function G has the form 

G{x,x'') = — \og\x-x^\+G{x,x''). (4) 

The harmonic function G is determined by the requirement that there is no flow in or out 
of the system. This implies the Dirichlet type boundary condition G = const, along the 
boundary or boundaries of the system. 

After inserting (p and A into the hydrodynamical Lagrangian and forming the equations 
of motion from the usual action principle, we obtain the following Lagrangian describing the 
point vortex system 

^ ^vortex ^^int ^^pot^ (5) 

where the vortex system interaction L'°* and irrotational Lp°* parts are 

^vortex _ ly"q''p{x'')e'^X\x'' 

a 

a 

-^E^VpMG'(a.^a;^), (6) 

L'°* = 5^ / d^xq''Jf°^e'^djGix,x'') (7) 

LP°* = L[p,^], (8) 

where G{x°') = G{x°-, x"') and JP°* is the part of the current density associated with irrota- 
tional flows. The relation of the Green's function to the vortex dynamics obtained here is 
similar to the results of 



The vortices interact with each other through a potential function given by the appro- 
priate Green's function G. In ([8]) we have obtained the self-energy of a single vortex, which 
is the energy associated with the flows created in the system by the presence of the vortex, 
by means of a simple renormalization procedure. The renormalization is performed by re- 
moving the singular logarithmic part of the Green's function from when the indices 
a and h coincide. In this coreless treatment the usual ultraviolet divergence present when 
calculating the vortex energy by direct integration is taken care of by this renormalization 
procedure. We see that the kinetic energy (or self energy) of a vortex is proportional to 
G[x) evaluated at the position of the vortex. This is a direct concequence of the defining 
properties of the Green's function. 

The equations of motion for the vortices are simply 

p(x'^)i" = J(a;''), a = l,...iV. (9) 

The study of vortex motion thus reduces to the hydrodynamical problem of studying the 
flows in the system. When non-vorticical flows are absent, the problem is to find the suitable 
Green's function under the Dirichlet boundary conditions defined above and then solve a 
group of nonlinear differential equations defined by ^ above. 



III. VORTEX DYNAMICS IN THE ANNULUS 
A. Vortex Hamiltonians in the annulus 



From this point on in the current paper, we will assume that the particle density p 
takes a constant value po everywhere. It should be noted that this fixing of p along with 
the requirement of energy conservation means that there is no energy transfer between the 
vortex system and the rest of the fiuid. This also means that vortex-antivortex annihilation 
and other processes that affect the vortex configuration are not possible in our model. The 
vortices initially placed in the system will exist forever. 

The form for the Green's function in the annulus can be found eg. in jl^ . It can be 
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expressed as 



1 



1 log ri log r2 



2-71 log 7] 



M^^oginr2)/2 + {e,-e2)/2,r]) 



(10) 



where the -^i are the Jacobi theta fuctions. Here we have chosen the inner and outer radius 
of the annulus to be -Ri = r/^^^ and R2 = respectively, so that the ratio of the inner 

radius to the outer radius is given by the parameter < r/ < 1 appearing in the theta 
functions. There is no fear of confusion since the -RjS can be removed and reinserted easily 
using their relation to rj and dimensional analysis. An important quantity is the geometric 
mean of the inner and outer circle bounding the annulus. It has the value 



(11) 



when expressed in terms of rj defined above. 

Since we are interested in the dynamics of one- and two-vortex systems, we write the 
two- vortex Hamiltonian in a form that is easily derived from the above Green's function 



H 



2 ^ m2po 

2tx 



I I logfrir?^/^) log(r2r7-'^/^) logr? 

- log |ri - r2| + — — 

logr/ 2 



n=l 



1 „/r" 



Vir2)" \r" r2 



a=l,2 



(12) 



where the H]^ are the radially symmetric renormalized single-vortex Hamiltonians given by 



rji qIpo \ \og^ [rr]^'^) ^ log 



log?7 



boo „ ^ 00 n„ 

2 ^'^(1-^^'')^ f'^''' ^^n{l-rj 



(13) 



We will use the above Hamiltonians in this and the next section to study one- and two- vortex 
systems in the annulus. 



B. Single vortex in an annulus 



The behaviour of a single vortex is obtained simply by studying the single- vortex Hamil- 
tonian (fT3|l . The result is a vortex moving with the renormalized flows it has created in the 
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Figure 1: Figure (a) shows R2 and between them the stationarity radius i?s as a function of 
q = R1/R2, all approaching the dashed line r = 1 as g — > 1. Figure (b) shows a typical graph of 
the single-vortex energy as a function of the distance from the center of the annulus. 



annulus 



V 



self ? rlog(rr/i/2^ 



2tx I r log 7] 

1 /(rr/V2)4n_^2n 
n=l \ I ) 

This implies the well known fact that a single vortex in an annulus will move along a circle 
at a constant distance from the origin jisl. We see that there is always a single value of r 
such that the vortex will remain at rest on the circle defined by that value. We will call this 
value of r the (single- vortex) stationarity radius and denote it with i?s. It is easily seen that 
i?s > 1 for all values of rj. The approximate value of i?s can easily be solved numerically 
(P (b)) for different rj. From the form of the single- vortex Hamiltonian ( |T3l ) we see that the 
circle r = Rs corresponds to a degenerate minimum value of the vortex energy ([D (b)). 

A single vortex cannot exist in the annulus if there is no angular momentum Lz in 
the system. We must thus also take into account the energy associated with this angular 
momentum. It can easily be shown that the vortex configuration that minimizes the energy 
of the fiuid, is the one that minimizes the free energy defined by 

F = E-ljLz, (15) 

where u is the angular velocity, with which the annulus is being rotated [l^. Since we are 
dealing with the system in a classical manner, the moment of inertia has the (classical) rigid 
value 



where N is the total number of particles in the gas 



ISj. The angular momentum of a vortex 



at distance r from the origin is easily calculated to be 



2 1] 4 ?7log?7 
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We will not get into the details of the equation for the minimum of the free energy. It suffices 
to note that, since cuL^ is a decreasing function of r, the minimum of free energy is located 
between Rs and the outer radius of the annulus. 



IV. VORTEX PAIRS IN THE ANNULUS 



We next proceed to study the dynamics of two- vortex systems in the annulus analytically 
and numerically. In subsections A and B we study analytically vortex-vortex and vortex- 
antivortex pairs respectively. In subsection C we study the motion of vortex pairs in the 
annulus numerically and demonstrate the results of subsections A and B. 



A. A vortex-vortex pair in an annulus 

It is easy to see that when the condition 9i — 62 = n is met, the motion of the vortex- 
vortex pair has simple analytic solutions for which 9i — 62 — 0. These type of solutions are 
obtained choosing the initial positions for the vortices as ^i(O) — ^2(0) = tt and ri(0) = '"2(0). 
This corresponds to a vortex-vortex pair moving on the same circle with constant angular 
velocity. The angular velocity of the corotating pair is 

1 log(r) 



/icorot ^ 



log(?7) 

n=l ^ 



1 1 

4n I ^4 ^2 



(18) 



Setting Og°=l^ = we see that there always exists a single solution for which the vortex pair 
remains at rest. This is analogous to the case of the single vortex. However, this time the 
stationary solution appears for a value r < 1 for the radii of the vortices. We will call this 
radius the vortex- vortex stationarity radius and denote it with i^g"*" (or -Rg", depending 
on the sign of the winding numbers). We have thus shown that for a vortex- vortex pair 
in an annulus there always exists infinitely many corotating solutions and also a stationary 
solution. The stationary solution minimizes the energy of the vortex pair. 

A solution satisfying the different condition ^1 -|- ^2 = is called counterrotating as the 
angular velocities of the vortices have opposite sign in this case. Such solution can appear 
when the vortices are on different sides of the stationarity radius; the stationarity radius 
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corresponds to a solution which is both corotating and counterrotating. The counterrotating 
solutions are not as easily described analytically as the corotating solutions, but they are 
easily studied numerically. 

As in the case of a single vortex, a vortex- vortex pair cannot exist in the annulus without 
the precense of angular momentum. The free energy is easily calculated since is linear 
with respect to the velocity fields and thus additive. Again, it suffices to note that the 
most favourable solutions are corotating solutions where both vortices are displaced (at 
least approximately) by the same amount from the stationarity circle towards the outer 
radius given the symmetry of the system and the complex nature of the vortex trajectories 
in solutions where ri and r2 differ by a "larger" amount (we will not study the concept of 
"large" in this context here). This is in accordance with the results of [20!| for the stable 
vortex configurations in an annular BEC. 



B. A vortex-antivortex pair in an annulus 

Because of the change of sign in the interaction part of the Hamiltonian, the conditions 
Oi — 62 = 'K and ri = r2 describe a counterrotating solution. Conversely, the solutions 
describing corotating behaviour are not easily expressed analytically but are easily studied 
numerically. Searching for a stationary solution leads us to study the equation 

i-4y ^ , , =0. (19) 

n=l ' 

This equation always has a solution r < 1 and thus there always exists a stationary vortex- 
antivortex -solution, this solution being again the solution which is both corotating and 
counterrotating and correspomding to a local minimum of the vortex-antivortex energy. We 
will call this value of r the vortex-antivortex stationarity radius and denote it with -Rg ~. 

It may be noted that for the case of the vortex-antivortex pair, there need not be any 
angular momentum in the system. In this case the free energy is just the energy of the 
vortex pair. One way to produce a vortex-antivortex pair confined in an annular region is 
to produce a soliton of topological charge in a BEC confined in a "mexican hat"-potential. 
This soliton can then decay, typically producing a vortex-antivortex pair on opposite sides 
of the annulus. Of course, such a soliton will carry with it some angular momentum. The 
situation still remains different from the case of one or two vortices of winding numbers of 
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the same sign. The qualitative properties of the corotating and counterrotating solutions 
suggest that they both can be at least metastable energetically. A confirmation for this fact 
can be obtained through direct numerical simulations of a ring-shaped BEC starting from 
the Gross-Pitayevskii equation. 



C. Numerical simulations 



The equations of motion describing vortex pairs in the annulus can easily be solved 
numerically using a standard adaptive stepsize Runge-Kutta method. We will choose the 
value of ?7 to be 1/3, since then the annulus is similar to that used in [13| (scaled by a factor 
of ~ 1.1547), where the inner radius was chosen to be 0.5cm and the outer radius 1.5cm. 
We notice that very important quantities, as expected, are the two- vortex stationarity radii 
i?g^ < 1 and i?g ~ < 1, which in our case have the numerical values R^^ ~ 0.74911 and 
^ 0.888928. 

Figures [2] (a-e) are the same cases as in [l^ in the figures 1-5, respectively. Figures [2] 
(e-f) show how sensitive the solutions are to initial conditions. In figures [3] and [4] we see the 
effect of the stationarity radii for the two- vortex solutions. The dashed circle in these figures 
is the circle r = -Rg for a vortex- vortex pair and the circle r = for a vortex-antivortex 
pair. 



D. A brief heuristic discussion 



The fact that we have been able to find analytic and numerical solutions with similar 
properties for both kinds of vortex pairs is easily explained by the geometry of the annulus. 
The annulus is not simply connected and the annulus has two boundaries. The first property 
has a profound consequence most simply illustrated by the vortex-vortex pair. We may 
think that the interaction between vortices is transmitted by the fiows of the original system 
(BEC, superfiuid or something else), or more simply that the interaction is transmitted 
only where the system exists (p 7^ 0), in this case the annulus. Thus we may think that 
vortex 1 interacts with vortex 2 through two different routes, "striking" at vortex 2 from 
both directions 6 O2— and 6 — > 6*2+. With this thought in mind it is not surprising that 
vortices on exactly opposite sides of the origin corotate or stay still, the interactions from 
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9— and 9-\- essentially cancel each other and all that is left is a self-induced velocity for the 
vortex pair. 

The two boundaries of the annulus have the following consequence: The harmonic part 
of the Green's function must satisfy the symbolic relation G{x^, x^) = G{x^) = oo on both 
boundaries. This and the finiteness of G inside the annulus directly implies that G{x^) 
must have a minimum value that corresponds to a stationary vortex solution. Because of 
the radial symmetry of the annulus this solution is infinitely degenerate and, as seen before, 
corresponds to a single circle cosentric to the circles bounding the annulus. The existence of 
this circle leads to the fact that the sign of the angular velocity of a vortex around the origin 
depends on whether it is inside or outside the stationarity circle. Adding more vortices to 
the system does not change the fact that the single vortex energies become infinite at the 
boundaries and that this leads to the existence of a stationarity radii. Thus vortex-vortex 
-pairs and vortex-antivortex -pairs can display similar behaviour. 

V. CONCLUSIONS 

We have shown that taking into account the boundary conditions of the doubly connected 
annulus leads to the existence of stationarity radii which greatly affect vortex motion in the 
annulus. This leads to the existence of corotating and counterrotating as well as stationary 
solutions for both vortex-vortex and vortex-antivortex pairs. 
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boundary changed, (d) Another vortex-antivortex case with both vortices initially inside the circle 
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performs motion similar to the vortex-antivortex pair in figure [2] (d). (f) A vortex-antivortex pair 
with initial positions inside and outside the circle r = performs motion similar to the vortex- 
vortex pair in figure [2] (b). 

[10] D. M. F. Chapman, J. Math. Phys. 19, 1988 (1978). 

[11] D. Crowdy and J. Marshall, Proc. R. Soc. A 461, 2477-2501, (2005). 

[12] D. Crowdy and J. Marshall, Phys. Fluids 17, 056602, (2005). 

[13] T. I. Zueva, Low Temp. Phys. 26, 85 (2000). 

[14] R. Courant and D. Hilbert, Methoden Der Mathematischen Physik I (Verlag Von Julius 

Springer 1924) p. 312. 
[15] V. L. Berezinsky, Zh. Eksp. Teor. Fiz. 61, 1144 (1971). 



13 



a b 




Figure 4: In figures (a) and (b) we see the change in the behaviour of a vortex-vortex pair when 
the initial position of the inner vortex is moved from outside the circle r = R^'^ to inside it. In (c) 

and (d) we see the corresponding change for a vortex-antivortex pair and the circle r = ■ (e) 
A corotating vortex-antivortex pair, (f) A vortex-vortex pair with the same initial positions as the 
vortex-antivortex pair in (e). 
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